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Abstract —Fifth-generation cellular networks are expected to 
exhibit at least three primary physical-layer differences relative 
to fourth-generation ones: millimeter-wave propagation, antenna- 
array directionality, and densification of base stations. In this 
paper, the effects of these differences on the performance of 
single-carrier frequency-domain multiple-access uplink systems 
with frequency hopping are assessed. A new analysis, which is 
much more detailed than any other in the existing literature 
and accommodates actual base-station topologies, captures the 
primary features of uplink communications. Distance-dependent 
power-law, shadowing, and fading models based on millimeter- 
wave measurements are introduced. The beneficial effects of 
base-station densification, highly directional sectorization, and 
frequency hopping are illustrated. 

1. Introduction 

This paper provides a performance analysis of the uplink 
of a fifth-generation cellular network in its primary mode 
of operation. It is assumed that millimeter-wave frequencies 
will be adopted, but that the basic structure of the fourth- 
generation (4G) single-carrier frequency-domain multiple- 
access (SC-FDMA) uplink systems will be maintained. The 
first assumption is widely supported in the current literature 
on fifth-generation (5G) systems Q, (a. The second one is 
based on the favorable characteristics of SC-FDMA and the 
critical importance of transmitter power efficiency, which is 
enabled by the low peak-to-average power ratio of SC-FDMA 
and similar single-channel modulations ID- It is further 
assumed that 5G systems will exploit frequency hopping, as 
is done in 4G/LTE systems fSl- 

The analysis in this paper applies the methodology of 
im, which we call deterministic geometry. Unlike stochastic 
geometry, deterministic geometry can accommodate arbitrary 
topologies with distance-dependent propagation models. The 
distance dependence of the propagation and fading models 
accounts for the fact that mobiles close to the base station 
have a line-of-sight (LOS) path, but the more distant mobiles 
do not. 

The conditional outage probability of a SC-FDMA uplink is 
derived, where the conditioning is with respect to an arbitrary 
network topology. In the numerical examples, we use an actual 
deployment of the base stations. Each simulation trial gener¬ 
ates a network realization in which the mobile placements 
are drawn from the uniform clustering distribution with each 
mobile having an exclusion zone 06). For each realization of 
the network, the outage probability is computed for a reference 
link. By averaging over many network realizations, the average 
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Fig. 1. Actual base station locations from a current cellular deployment. 
Base stations are represented by large circles, and Voronoi cell boundaries 
are represented by thick lines. 

outage probability and other statistical performance measures 
are computed. 

The remainder of the paper is organized as follows. Section 
ini describes the system model, which accounts for distance de¬ 
pendent fading, path loss, and shadowing based on millimiter- 
wave measurements, and for highly directional antennas and 
sectorization. Section III provides a closed-form expression 
for the conditional outage probability. Section IV provides 
a numerical evaluation of a 5G system. Finally, the paper 
concludes in Section V. 

11. Network Model 

In the network model, C base stations and M mobiles are 
confined to a finite area. As an example. Fig. [T] depicts an 
actual deployment of C = 121 base stations extracted from 
a database of base-station locations in the United Kingdom. 
Base stations are represented by large circles, and Voronoi 
cell boundaries are represented by thick lines. The network 
occupies 900 km^ inside a square. 

For practical reasons such as the high propagation loss 
and the complexity of channel state estimation at millimeter- 
wave frequencies, the antenna arrays at the base stations and 
mobiles exploit beamforming rather than spatial multiplexing 
for uplink transmissions (T) . Densification, high mobility, and 















the severe impact of blockages at millimeter-wave frequencies 
cause frequent handoffs and hence the need for rapid beam 
alignments. Sectorization, which is the division of base-station 
coverage into ( fixed angular sector beams centered at the base 
station, is used to reduce beam-alignment delays and pilot 
contamination. At millimeter-wave frequencies, the beams 
can be implemented using many antenna elements, perhaps 
hundreds, and hence have narrow beam widths and very small 
sidelobes and backlobes. Each mobile transmits through an 
adaptive and highly directional antenna array. 

The scalar Si^ I = represents the Ith sector 

or its receiver, and the scalar i = 1,2, represents 

the ith mobile. The vector S/, / = 1,2,...,(^C, represents 
the location of the Ith sector receiver, and the vector X^, 
i = 1, 2,..., M, represents the location of the ith mobile. The 
normalized sector beam pattern associated with Si is 


f 1, <0 s- 27r/C 

I 6, otherwise 


( 1 ) 


where is the offset angle of the beam pattern, and b is the 
relative sidelobe and backlobe level. 

Let Xi denote the set of mobiles served by sector Si. Let 
Xr G Xj denote a reference mobile that transmits a desired 
signal to a reference receiver Sj. Let g(i) denote a function 
that returns the index of the sector serving Xi so that Xi G Xi 
if g (i) = /. The sector that serves mobile Xi is assumed 
to be the one with minimum local-mean path loss when the 
mainlobe of the transmit beam of Xi is aligned with the sector 
beam of . Thus, the serving sector has index 


g(i) = argmax; (||S; - Xi||), Xi e Ai^ (2) 

where ^i^i is a shadowing factor for the link from Xi to Su 
/ (•) is the area-mean path-loss function, H-H is the Euclidean 
norm, and Ai denote the set of mobiles covered by the sector 
beam of Si. In the absence of shadowing, the serving sector 
will be the receiver that is closest to Xi. In the presence of 
shadowing, a mobile may actually be associated with a sector 
that is more distant than the closest one if the shadowing 
conditions are sufficiently better. 

The area-mean path-loss function is a function of the 
distance d between a source and destination and is expressed 
as the attenuation power law 

/ d\ 

^ d> do (3) 

where a (d) is the attenuation power-law exponent, and do is a 
reference distance that is less than or equal to the minimum of 
the near-field radius and the exclusion-zone radius of a mobile. 

The distance-dependent model of the power-law exponent 
for millimeter-wave frequencies takes into account the area- 
mean attenuation due to blockages and refiections that occur 
over network links. This model refiects the empirical fact 
that a (d) differs substantially for LOS and non-line-of-sight 
(NLOS) links, tending toward amin for the usually shorter 
LOS links and tending toward a much larger amax for the 


usually longer NLOS links 10, 0, Col- Empirical data 
indicates that usually there is a small range of link lengths 
for which there are both LOS and NLOS links. Therefore, 
a (d) is modeled as a monotonically increasing function: 

a (d) = amin + (o^max “ Q^min) tauh (/id) (4) 

which indicates that amin < a (d) < amax- The parameter 
/i controls the transition rate from amin to a value close to 

O^max* 

Local-mean large-scale terrain effects are included in the 
shadowing model. The shadowing factor can be derived from 
a deterministic terrain model or can be modeled as a random 
variable in a statistical model. In this paper, we assume 
lognormal shadowing in which the are independent, 

identically distributed, zero-mean Gaussian random variables 
with a distance-dependent variance. 

Lor millimeter-wave frequencies, empirical data El, El, 
ifTOl indicates that the effect of the shadowing increases for 
the usually longer NLOS links. Since there is a small range 
of link lengths for which there are both LOS and NLOS links, 
the standard deviation of the shadowing factor for millimeter- 
wave frequencies is modeled as a monotonically increasing 
function: 


(^) — <^min (<^max <^min) tunh (/id) (5) 

which indicates that (Jmin < cts (d) < (Tmax- 

The fading is assumed to have a Nakagami distribution 
function. Since the fading becomes more severe for the longer 
links, the distance-dependent Nakagami parameter is modeled 
as a monotonically decreasing function: 

m (d) = TTimax “ (^rimax “ ^min) tauh (/id) (6) 

which indicates that rumin <m{d)< rnmax- 

Lrequency hopping ifTTl is used in SC-LDMA uplink sys¬ 
tems to provide the diversity that will mitigate the effects 
of frequency-selective fading and inter sector interference. Be¬ 
cause of network synchronization and similar propagation de¬ 
lays for the mobiles associated with a cell sector, synchronous 
orthogonal frequency-hopping patterns can be allocated so that 
at any given instant in time, there is no intrasector interference. 
The frequency-hopping patterns transmitted by mobiles in 
other sectors are not generally orthogonal to the patterns in 
a reference sector, and hence produce intersector interference. 
The varying propagation delays from the interfering mobiles 
cause their frequency-hopping signals to be asynchronous 
with respect to the desired signal. Duplexing prevents uplink 
interference from downlink transmissions. 

Each mobile uses a frequency-hopping pattern over a hopset 
with L disjoint frequency channels. Let L/, I = 1,2,..., (^C, 
denote positive integer divisors of L such that L/Li > 2. Each 
mobile in A) is assigned a distinct block of Li contiguous 
frequency channels during each of its hop intervals, and the 
block may change to any of LjLi disjoint spectral regions with 
every hop. Consider an uplink reference signal that traverses 
a reference link from a reference mobile Xr to a reference 
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Fig. 2. Illustration of co-channel interference at base station Sj produced 
by an interfering mobile Xi arriving at the base station with a relative timing 
offset of tij. Each block in the illustration represents a hop, and the selected 
channel is indicated by the shading of the block. A partial-slot collision is 
indicated, where the interfering mobile has selected the channels used by the 
reference mobile during the second slot of the subframe. However, due to the 
asynchronism, the collision only extends for part of the slot. 


receiver Sj. Because of a possible incomplete spectral overlap, 
the received interference power from mobile Xi at Sj when the 
mobile’s signal collides with the reference signal is reduced 
by the spectral factor 

F; =min(Lj/i,,l). (7) 


Each active mobile points its antenna beam toward its base 
station. The beam pattern is modeled with two gains: one for 
the mainlobe and another for the sidelobes and backlobes. The 
normalized beam gain at the reference receiver Sj due to the 
angular offset of this beam pattern, is defined as 




1 , 


||S,-X.||||Sg(i)-X.|| > cos 1^2 J 
otherwise 


( 8 ) 


where 0 is the beamwidth of the antenna mainlobe, superscript 
T denotes the transpose, and a is the sidelobe and backlobe 
level. 

Associated with each potentially interfering mobile is a 
hop transition time ti^j at Sj relative to the hop transition 
time of a pair of hop intervals of the reference signal. Fig. 
[2] illustrates the relative timing between the reference mobile 
and the interfering mobile. The reference mobile transmits 
a turbo codeword of duration 2T, which is aligned with 
the subframe. It is assumed that the frequency separation of 
the two frequency channels of the two slots is sufficient for 
independent fading of fixed amplitude in each slot. If the base 
stations and mobiles are synchronized, then 


ti^j = [(||S, -X,|| - ||Sy -X,||)/c]modT (9) 


where c is the speed of an electromagnetic wave. As illustrated 
in Fig. [2l the reference signal encounters four time periods of 
potential interference from an active mobile 0 < t < 
ti^j ^ t < T < t < tij + T, and tij S T < t < 2T. 
The generic index k G {1,2, 3,4}, denotes a time period of 
the subframe with duration that varies with each Uj. In the 
example provided in Fig. [2l co-channel interference occurs 


during the third time period, where the interfering mobile 
has selected the channels used by the reference mobile. The 
fractional duration of the each of the four subframe time 
periods relative to the subframe period 2T are 






2T ’ 
2T 


k = 1,3 
k = 2,4 


( 10 ) 


The set of indices of potentially interfering mobiles is 
S = {i \ Xi ^ Xj^. Let Xi denote the number of mobiles 
associated with sector Si. Because of the required orthogo¬ 
nality of frequency blocks assigned to mobiles within each 
sector, Xi < LjLi and any additional mobiles within the 
sector are reassigned to other sectors. In view of the po¬ 
tential spectral overlaps, the maximum number of interfering 
mobiles within a sector during a subframe time period is 
min[max (Ly/L/, 1), X/]. Let Sk C S denote the set of 
interfering mobiles during subframe time period k. If Xi < 
max(I/y/I//, 1), then all Xi mobiles in sector I are in Sk- If 
Xi > max(I/y/I//, 1), then some of the mobiles in sector I 
cannot cause interference during subframe time period k. In 
that case, we approximate by randomly selecting a subset of 
the Xi mobiles to be included in Sk- 

Let Qi^k denote the probability that the signal from a 
potentially interfering mobile collides with the reference signal 
during subframe time period k, 1 < k < 4. The activity 
probability pi is the probability that mobile Xi is transmitting 
during time interval [0,2T). Assuming uniformly distributed 
frequency-hopping patterns that are orthogonal within each 
sector. 


(li,k = 




Li 


L 


-Pi^ 


i G Sk, 1 < k < 4 (11) 


and gi,/c = 0, otherwise. 

The instantaneous signal-to-interference-and-noise ratio 
(SINR) at sector receiver Sj when the desired signal is from 
Xr G Xj fluctuates because potentially interfering signals 
do not always coincide with the reference signal in time or 
frequency. Pilot sequences are used to estimate the complex 
fading amplitudes in the receiver. Therefore, the performance 
of the reference receiver is primarily a function of the average 
SINR defined as the ratio of the average power of the signal 
to the average power of the noise and interference, where 
the average is over the time interval of a subframe and turbo 
codeword. Thus, the average SINR during a subframe is 


7r,j — 




S/c=l h,kPid,kCiJ,k 


( 12 ) 


where M is the noise power, f^ j is the average received power 
from reference mobile X^, and pij^k is the received power 
from an interference signal that collides with the reference 
signal during subframe time period k. The indicators f^k are 
Bernoulli random variables with probabilities 


P[f,k — 1 ] — Q.i,k’) P[f,k — 0 ] — 1 Q.i,k 
i e Sk, 1 < k < 4. 


( 13 ) 



















Let Qij^k denote the fading gain of the signal from mobile 
Xi at Sj during time interval k. Assuming that the bandwidths 
of the L/Li and LjLj disjoint spectral regions exceed the 
coherence bandwidth, the are independent for each 

hop with unit-mean, and Qi^j^k — where ai^j^k has 

a Nakagami distribution with distance-dependent parameter 
rriij. Let Pi denote the maximum power from Xi that could be 
received at the reference distance do in the absence of fading 
and shadowing. Allowing for the spectral and beam factors, 
the received power from Xi at Sj^i G Sk^ during time interval 
k is 


(l|Si - Xi 

i ^ Sk, I < k < 4 




(14) 


where Oij is the arrival angle at Sj of a signal from and 
Bj {Oij) is the gain of the uplink beam pattern of Sj. 

Let gr,j,i and gr,j ,2 denote the unit-mean power gains due 
to the independent fading of the frequency-hopping reference 
signal in subframe slots 1 and 2, respectively. The power gain 
of independent Nakagami fading with parameter mo = mrj in 
each slot has the gamma density function: 


fr,j{x) = 


^ exp {—mox) 
r (mo) 


,{x) 


(15) 


where u{x) = 1, x > 0, and u{x) = 0, otherwise. The 
average power gain due to fading is g^ j = {gr,j,i + dr,i, 2 ) /2. 
Using (O, we obtain the density function of g^ j : 

(2mo)^'"°a;2™“"^exp(-2moa;) 

which is the power of a Nakagami random variable with 
parameter 2mo. This doubling of the Nakagami parameter 
indicates the beneficial effect of the frequency hopping in 
mitigating frequency-selective fading. The average received 
power from reference mobile X^ is 




where dr = 11S j 11 is the length of the reference link. 

Power control ensures that the local-mean power received 
from each mobile in a sector is equal to a constant. If this 
constant is the same for all sectors, then 


(||Sg(,) - Xiii), 

i £ Sk- (18) 


Substituting (m, dnll, and dH into (US, we obtain 

_ 9 r,j 

Lq + S/c=i '^ieSk 


(19) 


where 


^i,j — 


10[^^-i-^dg(o]/io/(||S. -X^ 
/(llSg(i)-Xill) 


is the ratio of the interference power from Xi to the reference- 
signal power, and 


ro = ( 21 ) 

is the signal-to-noise ratio (SNR) at the sector receiver when 
the fading is absent. 

III. Outage Probability 

Let [3 denote the minimum average SINR required for 
reliable reception of a signal from Xr at its serving sector 
receiver Sj^ j = g(r). An outage occurs when the average 
SINR of a signal from Xr falls below (3. The value of /3 sets 
a limit on the code-rate R of the uplink, which is expressed 
in units of bits per channel use (bpcu), and depends on the 
modulation and coding schemes, and the overhead losses due 
to pilots, cyclic prefixes, and equalization methods. The exact 
dependence ofRon/3 can be determined empirically through 
tests or simulation. 

The set of {^ij}, i ^ Sk^ for reference receiver Sj is 
represented by the vector flj . Conditioning on flj , the outage 
probability of a desired signal from Xr G Xj that arrives at 
Sj is 

e = P[7r,i </3|f^i] ■ (22) 

Because it is conditioned on Ptj, the outage probability 
depends on the particular network realization, which has 
dynamics over timescales that are much slower than the fading 
or frequency hopping. We define 

/3o = 2/3mo, ^ = rp (23) 

where the Nakagami parameter mo = for the reference 

uplink signal is assumed to be a positive integer. A derivation 
similar to the one in 0 yields 

2 mo —1 s 

e=l-e-P°^ E (24) 

5=0 t=o 

where 

4 

E n n (25) 

£ik>0 k=lieSk 

Efc=i Eie5^ ^ik=t 

the summations in ESi are over all sets of nonnegative indices 
that sum to t, 




e = o 


Ge{t,j,k) - qi,krie+mij) (^nnj+e 

' V rriij J 


£\r{mij) 


KX ^>0 


(26) 


and 


+ i^Sk, l<k<4. (27) 


( 20 ) 










IV. Numerical Results 

In the following examples, performance metrics are calcu¬ 
lated by using a Monte Carlo approach with N simulation 
trials. In each simulation trial, a realization of the network of 
Fig. 1 is obtained by placing M mobiles within it according to 
a uniform clustering distribution with each mobile having an 
exclusion-zone radius set equal to do = 0.004 km. Randomly 
generated shadowing factors are used to associate mobiles with 
cell sectors and in other computations. The code rate permitted 
by the threshold is given by 

R = log2 (1 + lsl3) (28) 

where Is = 0.794 corresponds to a 1 dB loss relative to 
the Shannon bound for complex discrete-time AWGN chan¬ 
nels. The density of mobiles in a network of area Anet is 
A = M/Anet- The outage probability of the reference link 
for simulation trial i is computed by applying (l24l)— (|27]) . 
The throughput of the reference uplink for simulation trial 
i is R{1 — ei). The average outage probability over all the 
simulation trials is 

1 ^ 
i=l 

The maximum rate of successful data transmissions per unit 
area is characterized by the area spectral efficiency, defined 
as 

A=XR{l-e) (30) 

where the units are bits per channel use per unit area. 

To avoid edge effects, the performance is measured for 
mobile stations placed in the yellow shaded area of Fig. 1, 
which is a 20 km by 20 km square located in the mid¬ 
dle of the network. To consider the effect of base-station 
densification, the base-station deployment and the density of 
mobiles A = 20Ikmf are maintained while the size of the 
network is scaled (dimensions redefined). For a fixed number 
of base stations, each scaling reduces the associated number of 
mobiles and decreases the transmission distances accordingly. 
In considering the effects of intersector interference, only the 
strongest 30 signals were used, as the attenuation of further 
signals was severe enough to make their individual effects 
negligible. 

The slot duration is T = 0.5 ms, PrfAf = 30 dB, 
L/Lj = L/Li = 10, and p = 40/km. Each base sta¬ 
tion has b = 0.01 and ( = 24. Each mobile beam width 
is 0 = 0.1 TT radians and sidelobe level is a = 0.1, 
and the common activity factor is pi = 1. The propaga¬ 
tion parameters are (rTlrnax?U^min? U^max? <Tmax) — 
(2,1, 2.3,4.7,6.1 dB, 12.6 dB), which are suitable for an ur¬ 
ban network. The reference signal has fading parameter mo = 
1 when dr = 0.1 km, and mo = 2 when dr = 0.01 km. The 
ratio C/M serves as a measure of the densification of the base 
stations. 

To capture the impact of densification, we need to take 
into account the decrease in the typical link-length dr of the 
reference link and the consequent increase in Fq as C/M 



Eig. 3. Average outage probability for dro =0.1 km, /3 = 3 dB, N = 10^, 
and distance-dependent fading. 



Fig. 4. Area spectral efficiency for three values of the SINK threshold, 
dro = 0.1 km, and N = 10^. 

increases. For a fixed density of mobiles, the average cell 
radius and hence a typical dr are proportional to 1/ ^C/M. 
For the range of C/M of interest, we consider 

^- ^ ’-0 0.01 < C/M < 1 (31) 

where dro is the typical link length when C/M = 0.01. 

Fig. [3] depicts the average outage probability e for a typical 
link as a function of the densification for dro = 0.1 km, (3 = 3 
dB, N = 10^, distance-dependent fading, both shadowing 
and no shadowing (cTmin = cTmax = 0), and both frequency 
hopping and its absence. In the absence of frequency hopping, 
the reference signal is assumed to experience constant fading 
over a subframe, which is valid if subframe duration is less 
than the coherence time. The figure illustrates the reduced 
















Fig. 5. Outage probability for eight uplinks and the average outage 
probability over all the uplinks for CjM = 0.1, jS = 3 dB, distance- 
dependent fading, and a single simulation trial. 

importance of frequency hopping in the presence of substantial 
shadowing because of the decrease in the typical values of 
Fq. The primary reasons for the monotonic decrease in e with 
densification and the increased importance of the frequency 
hopping are the increases in Fq and the reduced shadowing 
and fading experienced by the reference signal. Calculations 
show that decreasing L/Lj has a very small detrimental 
effect, which indicates the minor importance of the intercell 
interference due to the sectorization and the narrow adaptive 
beams. 

Fig. [4| depicts the area spectral efficiency ^ as a function 
of the densification for three values of the SINR threshold, 
dro = 0.1 km, and N = 10^. Results are shown for both 
the distance-dependent fading and the more severe Rayleigh 
fading for which rrii^j = m{d) = mo = 1. Increases in the 
SINR threshold p of the network links increase the outage 
probability. However, for sufficiently large values of C/M and 
hence Fq, this effect is minor compared with increased code 
rate that can be accommodated. As a result, the area spectral 
efficiency increases significantly. 

Fig.ia shows the outage probability as a function of code¬ 
rate R iox C/M = 0.1, /3 = 3dB, distance-dependent 
fading, and a single simulation trial. The dashed lines in the 
figure were generated by selecting eight random uplinks in the 
network and computing each outage probability. The average 
over all the uplinks is represented by the solid line. Despite 


the use of power control, there is considerable variability in 
the dependence of the outage probability on the code rate due 
to the irregular network topology. 

V. Conclusions 

This paper derives an analytical model for calculating the 
outage probability and area spectral efficiency. The model in¬ 
cludes the effects of millimeter-wave propagation, directional 
beams, frequency hopping, an arbitrary network topology, and 
the assignment of frequency blocks to mobiles. Numerical ex¬ 
amples illustrate the effects of various features and parameters. 
Base-station densification improves the network performance 
significantly. The significance of the intercell interference is 
greatly reduced because of the highly directional sectorization 
and beamforming. The frequency hopping more effectively 
compensates for frequency-selective fading as the shadowing 
decreases. The area spectral efficiency increases with the code 
rate despite the increase in the outage probability. 
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